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HYPER  VELOCITY-IMPACT  SHOCK  INDUCED  DAMAGE  TO  STEEL  ARMOR 


I.  INTRODUCTION 

An  In-housc  Laboratory  Independent  Research  (ILIR)  program  was  conducted  at 
the  Mobility  Equipment  Research  and  Development  Command  (MERAl)COM)  to 
design  hypervelocity  projectiles  that  maximize  backfacc  spallation  in  steel  armor.  The 
projectiles  are  small  (0.5  to  5 g eaeh)  so  that  several  hundred  ean  be  located  on  one  end 
of  an  explosive  driver  weighing  15  to  30  kg.  The  projectiles  would  be  ejected  in  a near- 
vertical downward  direction  toward  a target  in  a shotgun  pattern  at  eject  velocities  of 
8 to  10  kin/s  and  at  a stand-off  distance  of  20  to  50  m.  This  enables  the  consideration 
of  large  lethal,  or  damage,  radii  of  the  order  of  tens  of  meters.  The  primary  target  for 
•the  experimental  program  has  been  steel  armor  2 to  5 cm  thick.  Damage  would  be  by 
backfacc  spallation  and  fragmentation  and  the  resultant  effects  on  personnel,  muni- 
tions, fuel,  equipment,  instruments,  and  soft  interior  components  (of  military  vehicles) 
such  as  cabling. 

When  the  program  started  two  years  ago,  there  was  insufficient  information  on 
spallation  physics  to  design  the  type  of  projectile  desired.  Accordingly,  a program  was 
startl'd  to  obtain  the  needed  experimental  and  analytical  information.  The  specific 
goal  was  to  design  a projectile  capable  of  shaping  the  pulse  induced  in  the  steel  target 
in  such  a way  that  a large  amount  of  steel  would  be  ejected  from  the  backfacc  of  the 
specimen  in  the  form  of  many  small  damaging  fragments.  During  the  course  of  the 
research,  new  phenomena  were  observed  and  arc  reported  herein.  In  particular,  the 
importance  of  the  130-kbar  phase  change  in  martensitic  steel  and  the  resultant  effect 
on  spallation  arc  discussed. 

Over  100  successful  impact  experiments  wen'  performed  in  the  Naval  Research 
Laboratory  (MIL)  light  gas-gun  facility.  The  experiments  wi  re  designed  and  the 
results  were  analyzed  at  MERADCOM.  Metallurgical  examination  of  the  impacted 
specimens  was  done  at  the  Stanford  Research  Institute  (SRI),  with  some  assistance 
from  Sandia  Laboratories  Albuquerque  (SLA).  The  2-d  code  calculations  were  done 
by  SLA.  Mechanical  material  properties  of  the  target  materials  were  measured  at 
MERADCOM. 

Ill  order  to  develop  reliable  theory  and  associated  analytical  tools,  simple  spheri- 
cal projectiles,  mostly  weighing  less  than  1 g,  were  used.  Impact  angles  were  generally 
normal,  although  a few  non-normal  impacts  were  made.  Also,  a few  impacts  were 
made  with  both  simple  and  composite  non-spherical  projectiles.  Projectile  materials 
included  nylon,  lexan,  water,  heavy  oil,  steel,  and  several  different  ceramics. 
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These  choices  of  projectile  materials  were  based  on  the  fact  that  the'  proposed 
warhead  is  weight-limited  rather  than  volume-limited.  This  enables  many  projectiles  of 
the  same  size  to  be  carried,  or  larger  projectiles  may  be  used.  Size  is  important  because; 
tin;  projectile  diameter  should  be  of  the  order  of,  or  greater  than,  the  target  thickness 
in  order  to  optimize  spall.1 2 *  The  vaporization  or  sublimation  energy  of  the  projectile 
material  should  he  small  compared  to  the  kinetic  energy  of  impact.  The  general  idea  is 
to  convert  most  of  the  impact  energy  into  induced  shock.  This  means  that  inefficient 
processes  such  as  penetration  should  be  minimized.  Most  of  the  targets  were  wrought, 
homogeneous  steel  plate  stock  1.25  cm  and  2.5  cm  thick.  The  steed  was  a carbon 
mangaiiese-iiickel-chromium-molybdenum  alloy  which  was  austenitized  at  900°  C, 
water  quenched,  and  tempered  at  540°  C.  The  representative  hardness  was  360  Urincll 
(Him),  and  the  approximate;  quasi-static  tensile  strength  was  12  kbars.  The;  dynamic 
spall  strength  corresponded  to  a critical  tensile  pressure  of  38  kbar.2  3 4 5 

After  impac  t,  the  target  specimens  were  dissected  through  the  crater  center, 
polished,  etched,  and  examined  metallurgically  by  SRI,  SLA,  and  MERADCOM. 
Extensive  2-d  code  calculations  of  the  impact  spall  phenomena  for  oik;  specimen  were4 
made  by  SLA.  SLA  id  so  measured  spall  strength  and  equation  of  state. 

The  experimental  results  are  summarized  in  Section  II.  They  art;  analyzed  and 
discussed  in  Section  III  The  2-d  computations  performed  by  SLA  arc  described  in 
Section  IV.  The  130-kbar  phase  change  and  spall  phenomena  art'  discussed  in  the 
remaining  sections. 


II.  EXPERIMENTAL  RESULTS 

A photograph  of  a dissected  impact  crater  is  shown  in  Figure  la.  The  target 
material  was  1.25-cm- thick  wrought  steel.  The  projectile  was  a 0.52-g  nylon  sphere, 
with  a diameter  of  0.95  cm  and  an  impact  velocity  of  5.18  km/s.  The  crater  is  repre- 
sentative of  craters  formed  by  impact  of  lexan-encased  I120  and  solid  nylon  spheres. 
The  walls  of  such  craters  exhibit  a series  of  evenly  spaced  ridges,  or  serrations,  which 
are  symmetric  about  a central  axis  of  rotation.  The  macrooracks  observed  extend 
downward  from  the  bottom  of  the  valleys  in  the  crater  floor.  In  the  case  of  steel-on- 
steel  impacts  (Figure  lb),  the  serrations  in  the  craters  are  not  observed  and  the  macro- 


^ There  is  some  experimental  evidence  that  much  smaller  (non -penetrating)  projectiles  are  efficient  spoliators.  (A 
spoliator  is  a projectile  that  optimizes  backfacc  spall.  A pcnctrator  optimizes  penetration.) 

2 

J.  W.  Bond  and  G.  W.  Ullrich,  “Two-Dimensional  Spallation  induced  by  Hypcrveloeily  Impact  in  Wrought  Steel 
Plate,”  USAMERDC  Report  1067 July  1973. 

^ A.  Stevens,  Personal  Communication,  SLA,  1973. 

* D.  A.  Shock cy,  “Support  of  Armor  Fragmentation  Studies,”  SRI  Project  PYU  215.1,  15  October  1973. 

° I).  A.  Shockey,  “Hypcrveloeily  Impact  and  Associated  Phenomena,”  SRI  Project  PYl)  2151,  18  December  1973. 
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Figure  2.  Phase  change  and  associated  damage  below  the  crater  resulting  from  the  impact  of  an  H2  O/Lexan  sphere  on  2.5-cm 
wrought  steel.  The  projectile  mass  was  0.32  g and  the  impact  velocity  was  6.03  km/s. 


eraeks  extend  upward.  (This  difference  has  not  been  explained.)  Figure  2 shows  a 
10X  photograph  of  the  bottom  of  the  center  of  a crater  produced  by  a 6.03-km/see 
impact  of  a 0.35-g  H20/Lexan  sphere  on  a 2.5-cm  wrought  steel  target.  The  apparent 
shaded  region  extends  completely  around  the  crater  and  is  readily  visible  on  almost  all 
of  the  impact  specimens.  This  region  has  extremely  high  hardness  (500  Bhu)  and  dis- 
plays a fine,  untempered,  martensitic  mierostrueture.  This  observed  metallurgical 
transformation  has  been  shown  to  be  due  to  the  130-kbar  polymorphic  phase  transi- 
tion in  iron  or  martensitic  steel.6  7 The  fine  adiabatic  shear  lines  are  readily  observable 
on  this  photograph;  also,  note  the  large  voids  below  the  bottom  of  the  crater. 

When  the  stress  wave,  induced  at  impact,  travels  through  the  target  and  reflects  in 
tension  from  the  baekfaee,  a spall  layer,  or  gap,  is  produced  (as  seen  in  Figure  1).  Spall 
is  not  observed  for  these  conditions  at  velocities  below  about  3 kin/s.  At  slightly 
greater  velocities,  incipient  spall  is  observed  by  the  appearanee  of  mieroeracks.  The 
microcracks  coalesce  to  form  a complete  spall  layer  at  about  3.1  km/s.  As  the  impact 
velocity  is  increased,  the  width  and  diameter  of  the  spall  layer  also  increase,  and  at 
some  higher  velocity  the  backfacc  plugs  or  fragments.  Non-adiabatic  shear  occurs  at 
the  outer  edge  of  the  spall  layer.  For  a 0.52-g  nylon  sphere  on  1.25-cm  wrought  steel, 
backfacc  plugging  occurs  at  about  5.8  km/s.  For  smaller  (lighter)  spheres,  or  for 
thicker  targets,  the  plugging  (or  baekfaee  fragmenting)  velocity  is  higher. 

In  spite  of  the  extensive  experimentation  just  described  (as  well  as  all  other 
experimentation  performed  in  the  U.S.),  it  is  still  not  possible  to  design  an  effective 
spallator.  A few  preliminary  guidelines  have  been  obtained,  however.  These  are 
complicated,  henee  there  is  insufficient  space  in  this  report  to  present  details.  Briefly, 
it  is  clear  that  pulse  shaping  is  the  key  to  the  production  of  damaging  baekfaee  spall. 
This  is  exemplified  in  Figure  3 which  shows  the  reflection  and  resultant  spall  of  a tri- 
angular stress  pulse.  The  incident  pulse  reflects  in  tension  at  the  free  surface  x0. 
When  the  spall  threshold  is  exceeded,  a spall  layer  is  formed  at  x^  The  distance  Xj- 
x0  must  be  sufficiently  small  so  that  the  backfacc  can  be  spalled  and  fragmented  but 
sufficiently  large  so  that  a damaging  amount  of  baekfaee  material  is  ejected. 

After  the  spall  layer  is  formed  at  xl9  the  subsequent  stress  pulse  reflects  at  x, 
rather  than  at  x0.  This  can  result  in  the  formation  of  a second  spall  layer  at  x2.  Ob- 
viously, this  sequence  of  events  can  be  repeated  to  form  multiple  spall  layers.  In  order 
to  prove  this  hypothesis,  several  different  projectile  configurations  were  tested.  One  of 
these  configurations  was  a hollow  plastic  sphere.  The  reasoning  here  was  that  upon 
impact  more  mass  would  be  delivered  at  the  impact  periphery  and  "double -impact” 
conditions  for  thin  flyers  would  exist  at  the  impact  center.  This  heuristic  prediction 

6 D.  A.  Shoekey,  “Support  of  Armor  Fragmentation  Studies,”  SKI  Project  PYU  2J51,  15  October  1973. 

n 

' D.  A.  Shoekey,  “Hypcrvelocity  Impact  and  Associated  Phenomena.”  SRI  Project  PYU  2151,  18  December  1973. 
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Figure  4.  Spallation  damage  in  1.25-cm  wrought  steel  by  hollow  Lexan  sphere.  Mass  = 
0.565  g;  diameter  = 1.08  g;  impact  velocity  = 5.6  km/s. 
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was  eminently  borne  out  in  several  experiments  in  which  the  projectiles  were  hollow 
spheres.  An  example  is  shown  in  Figure  4.  In  this  case  the  projectile  was  a hollow 
lexan  sphere  weighing  0.56  g with  a diameter  of  1.08  cm  and  an  impact  velocity  of 
5.59  km/s.  Note  the  double  spall  layer  and  the  backface  fragmentation. 

A few  other  impact  experiments  were  performed  in  which  the  projectile  was 
computer-designed  to  optimize  backface  spall.  These  experiments  were  particularly 
successful  in  showing  that  backface  spall  can  be  optimized. 

III.  DISCUSSION  OF  EXPERIMENTAL  RESULTS 

In  general,  a well-conceived  experimental  program  produces  results  that  have 
not  been  predicted  theoretically.  This  program  is  no  exception.  The  experimental 
“surprises”  are  too  numerous  to  cover  in  this  report;  but  a few  of  the  more  important 
results  can  be  discussed: 

1.  Serrations  on  the  crater  floor  for  nylon-  and  water-steel  impacts  have  not 
been  explained.  In  particular,  the  reason  for  the  serration  formation  on  nylon-steel 
impact  and  the  lack  of  serration  formation  on  stecl-stcel  impact  have  not  been 
explained. 

2.  Preliminary  or  tentative  explanation  of  the  iriacrocracks  observed  for  nylon- 
on -steel  impacts  has  been  given.  Hut  this  explanation  fails  to  show  why  the  macro- 
cracks extend  upward  for  steel-steel  impacts  as  compared  to  downward  for  nylon-on- 
stcel  impacts. 

3.  Shear  lines  are  observed  below  the  crater  floor  in  almost  all  of  the  cases. 
Their  role  in  the  damage  process  has  not  been  explained,  however. 

4.  The  macroscopic  voids  observed  below  the  crater  center  have  not  been 
explained. 


5.  The  importance  of  the  observed  130-kbar  phase  change  has  been  explained 
for  one  case.  This  is  discussed  in  Section  IV.  Much  more  work  is  needed  on  this 
phenomenon. 

6.  Spallation:  the  2-d  computer  calculations  performed  by  SLA  have  shown 

excellent  agreement  with  one  impact  experiment.  This  applies  to  the  configuration 
and  location  of  the  spall  layer  for  a relatively  simple  ease:.  If  this  work  could  be  ex- 
tended to  more  important  and  relevant  cases  it  would  be  of  extreme  significance. 
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7.  General  Summary  . 


a.  Metallographic/Fractograpliic  Features.  As  indicated  previously,  con- 
siderable circumstantial  evidence  leads  to  the  conclusion  that  the  dark-etching  hemi- 
spherical volume  of  material  under  the  impact  site,  such  as  shown  in  Figure  2,  has 
undergone  the  pressure -induced  phase  transformation: 

a(bcc)^e<I,cp) 

Material  within  this  zone  is  of  significantly  finer  grain  structure  and  higher  hardness 
than  material  elsewhere.  Back  surface  laser  interferometer  records  obtained  in  plate 
slap  experiments  by  SLA  exhibit  clearly  the  disturbance  attributable  to  a wave  re- 
flected from  a denser  (e)  phase.8 

Profuse  shear  bonding  was  observed  in  the  subcrater  region  and  is  felt 
to  play  an  important  role  in  the  cratering  process.  The  long  cracks  extending  into  the 
armor  steel  from  the  crater  walls  tend  to  lie  on  shear  bands  and  tend  to  join  with  iso- 
lated pieces  of  material  that  form  part  of  the  ejecta. 

Shear  bands  that  are  not  associated  with  the  crater  walls  are  also  promi- 
nent in  the  transformed  region.  A number  of  equiaxed  micro  fractures  having  a ductile 
appearance  are  generally  present  in  this  region,  linked  together  by  shear  bands,  which 
are  easily  observable  because  of  their  distinct  etching  behavior.  Again,  cracks  are 
observed  to  have  formed  in  the  bands  and  suggest  the  mechanism  of  material  removal 
in  this  region.  No  evidence  of  shear  was  found  in  the  back  surface  regions. 

Fracture  damage  in  the  back  surface  region  proceeds  by  the  nueleation, 
growth,  and  coalescence  of  microfractures,  followed  by  a widening  of  the  resulting 
liiaeroerack  and  subsequent  scabbing  of  the  back  surface  by  a shearing  process  that 
does  not  appear  to  be  adiabatic.  Observations  of  back  surface  damage  at  all  stages  of 
development  were  possible  because  of  the  wide  variations  in  impaet  conditions. 
Detailed  stress  histories  were  not  obtained,  however. 

b.  Correlations  with  Stress  History.  An  attempt  was  made  to  correlate 
inetallographic  and  fractographie  features  observed  in  the  specimens  with  stress 
histories.  Only  one  impaet  experiment  (described  in  Section  IV  of  this  report)  was 
simulated  computationally  (by  SLA).  Many  of  the  essential  features  of  the  stress  his- 
tory, however,  may  be  similar  in  other  experiments.  Thus  the  following  correlations 
may  be  common  to  hyperveloeity  impaet  in  general  (although  many  more  experiments, 
metallurgical  analyses,  and  theoretical  extensions  are  needed  before  this  statement  can 
be  properly  delineated).  The  observed  features  for  the  single  ease  described  here  were 
predicted  (post-test)  with  outstanding  accuracy  by  the  SLA. 

® I).  A.  Shockey,  “Hypervclocily  Impact  and  Associated  Phenomena,”  SKI  Project  PYU  2151,  18  December  1973. 
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(1)  The  boundary  of  the  distinct  etching  region  that  had  undergone 
the  polymorphic  phase  transformation  corresponds  to  a 130-kbar  isobar.  Furthermore, 
SLA  showed  that  the  effect  of  the  phase  transformation  is  to  increase  the  stress  pulse 
duration  and,  hence,  the  back  surface  damage.  Thus,  all  other  things  being  equal,  an 
armor  steel  that  does  not  exhibit  this  phase  transformation  (such  as  fee  steel)  should 
be  less  subject  to  impact  spallation  damage  than  one  that  does. 

(2)  The  ductile  void-like  microfractures  in  the  suberater  region  corre- 
late qualitatively  with  the  high  triaxial  tensile  stresses  and  the  elevated  temperatures 
that  exist  there  shortly  after  impact.  Shock  heating  occurs  as  the  compressive  wave 
propagates  through  this  region,  driving  the  temperature  up  and  relaxing  the  yield  stress. 
Release  waves  run  in  from  that  part  of  the  impact  surface,  not  under  the  projectile,  and 
intersect  under  the  crater,  creating  a stress  state  having  a large  component  of  spherical 
tension.  Mierofraetures  begin  to  nucleate  and  grow,  and  the  elevated  temperatures 
favor  the  ductile  void  growth  mode. 

(3)  The  mierohardness  profiles  are  also  in  accord  with  the  computed 
stress  histories.  A test  profile  made  on  the  section  surface  in  the  direction  of  impact 
showed  a high,  reasonably  constant  level  of  hardness  in  the  transformed  zone  that 
fell  sharply  to  the  original  hardness  at  the  boundary  of  the  dark,  or  shaded,  region. 
Profiles  taken  at  45°  to  the  impact  direction,  however,  exhibited  a double  humped 
behavior.  The  initial  drop  in  hardness  (within  the  transformed  zone)  detected  by  45° 
profiles  is  attributed  to  release  waves  that  run  from  the  free  impact  surface  into  the 
suberater  regions  and  partially  unload  the  material.  The  single  humped  profile  in  the 
impact  direction  is  consistent  with  this  being  a line  of  symmetry.9 

(4)  Shear  banding  in  the  suberater  region  is  also  consistent  with  the 
computed  stress  history,  high  shear  stresses  being  predieted  in  the  armor  near  the 
projectile- target  interface.  The  shear  bands  connecting  the  suberater  voids  are  thought 
to  form  rather  late  in  the  event.  The  original  spherical  tensile  field  is  distorted  by  the 
nucleation  and  growth  of  voids,  which  give  rise  to  significant  shear  stresses.  The  shear 
stresses  are  relieved  by  plastic  flow,  which  is  concentrated  in  narrow  bands  running 
between  the  voids.  The  material  within  the  bands  gets  very  hot  but  is  quenched 
immediately  by  the  steel  on  either  side.  Therefore,  this  material  transforms  to  aus- 
tenite (fee)  and  then  quickly  to  martensite,  which  is  hard  and  brittle,  and  acquires 
cracks  under  subsequent  stress  reverberations. 

(5)  The  formation  of  spallation  or  fraeture  damage  is  perhaps  the  best 
understood  feature  observed  in  the  impact  specimens,  as  evidenced  by  the  excellent 
agreement  with  the  SLA  predictions  of  damage  morphology.  The  damage  occurs  by 


^ O.  A.  Shockcy,  “Hypervelocity  Impact  and  Associated  Phenomena,'’  SKI  Project  PYIJ  2151,  18  December  1973. 
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the  nucleation,  growth,  ami  coalescence  of  microcraeks  when  the  compressive  pulse 
reflects  from  the  back  surface  and  interacts  with  the  release  wave  from  the  front  sur- 
face. The  formation  and  fragmentation  of  baekface  spall  is  not  well  understood, 
however. 


(6)  As  a result  of  the  effects  and  measurements  discussed  previously 
and  the  two-dimensional  calculations  discussed  in  Section  IV,  it  is  possible  that  the 
magnitude  and  duration  of  the  tensile  pulse  calculated  by  the  SLA  2-d  codes  can  be 
used  in  conjunction  with  the  material  dynamic  fracture  parameters  to  predict  the 
extent  of  back  face  fracture  damage  in  cjuantitative  microscopic  detail.  However,  more 
quantitative  correlations  clearly  require  more  experiments  and  more  detailed  2-d 
computations,  in  addition  to  quantitative  metallurgical  examination  of  the  specimens 
(counting,  measuring,  and  positioning  individual  artifacts  on  the  polished  surfaces). 

(7)  Measurements  show  that  at  equal  kinetic  energies,  spherical  nylon 
projectiles  are  as  effective  as  steel  projectiles  in  producing  back  surface  damage,  and  at 
equal  momentum  nylon  projectiles  are  much  more  effective  than  steel  projectiles.  It 
is  not  clear  at  this  time,  however,  which  is  the  more  important  parameter— energy  or 
momentum. 


IV.  2-d  COMPUTATIONS 

It  was  fortuitous  that  the  MERADCOM  impact  spallation  experiments  began 
about  one  year  before  SLA  scheduled  their  new  two-dimensional  hydrodynamic  codes 
for  usage.  Accordingly,  in  May  1973,  SLA  began  to  code  the  impact  of  a nylon  sphere 
on  1.25-cm  steel  armor.  The  experimental  results  are  showm  in  Figure  1.  The  impact 
velocity  was  5.18  km/s  and  the  projectile  mass  0.52  g. 

Two-dimensional  code  calculations  of  this  test  were  performed  at  SLA  using  both 
the  CSQ  Eulerian  Code  and  the  TOODY  Lagrangian  Code  with  rezoning.10  The  initial 
impact  pressure  was  about  400  kbar,  which  causes  iron  and  martensitic  steel  near  the 
impact  point  to  undergo  the  ot  e polymorphic  phase  change.  The  nylon  sphere 
greatly  distorts  and  goes  into  a liquid,  partially  vapor  state.  The  steel  plate  suffers  a 
large  distortion  near  the  impact  point  and  requires  both  an  accurate  elastic-plastic 
model  and  a good  material-failure  model  elsewhere  in  the  plate.  Also,  a high  degree  of 
resolution  is  required  for  predicting  complete  or  internal  spallation  because  the  atten- 
tuation  of  the  initial  compressive  pulse,  its  reflection,  and  the  interaction  of  the  com- 
pressive and  the  reflected  release  wave  must  all  be  calculated  accurately. 


10 


D.  A.  Shockcy,  “llypervelocity  Impact  and  Associated  Phenomena,’*  SRI  Project  PYU  2151,  18  December  1973. 
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The  130-kbar  phase  change,  the  correct  spall  strength,  and  an  adequate  failure 
model  have  all  been  determined  to  be  important  for  predictive  and  interpretive  calcula- 
tions. A detailed  examination  of  the  stress  histories  for  the  cases  with  and  without  the 
phase  change  has  shown  why  different  failures  should  he  expected  in  these  two  eases.11 
When  the  phase  change  is  included  in  the  calculations,  the  loading  portion  of  the  stress 
wave  consists  of  two  parts  as  the  wave  separates  at  130  kbar.  More  important  is  the 
difference  observed  as  unloading  occurs;  a rarefaction  shock  is  present  when  the  phase; 
change  is  included.  These  differences  in  die  stress  wave  structure  not  only  cause  a 
cylindrical-conical  failure  to  occur  directly  below  the  crater  but  also  to  result  in  a 
propagating  pulse  that  is  nearly  square  for  the  case  with  a phase  change;  and  nearly 
triangular  without.  Upon  reflection  from  a free  surface,  the  square  pulse;  transfers 
essentially  all  the  momentum  to  the  spall  layer,  whereas  the  triangular  pulse  is  not 
nearly  as  effective  in  momentum  transfer.  This  clearly  shews  the  importance  of  the; 
stress  pulse  shape;  in  producing  effective  spall. 

The  measured  spall  stress  of  38  khar12  13  for  the'  MS  12560  steel  was  used  in  the 
foregoing  calculations,  and  the  results  showed  excellent  agreement  with  the*  experi 
rnent.  This  was  especially  true  for  the  Lagrangian  results  where  the  crater  diameter, 
crater  depth,  spall  layer  thickness,  spall  length,  spall  bulge,  cylindrical-conical  fai lure, 
and  the  area  that  undergoes  a phase  change  all  show  nearly  one-to-one  correspondence 
with  the  experiments.  The  calculations  did  not  re  product;  the  serrations  in  the  crater 
floor,  however,  nor  the  maerocracks  and  shear  hands  observed  below  the  crater  floor. 


When  the  phase  change  is  included  in  both  tin;  CSQ  and  the  TOODY  calculations, 
the  only  significant  difference  is  the  spall  length  (or  diameter).  This  difference  is  a 
direct  result  of  the  material  failure  treatment  at  the  spall  plane  where  a material  coordi- 
nate treatment  allows  a higher  accuracy.  Tin;  very  good  numerical  and  experimental 
agreement  demonstrates  tin*  capability  of  the  codes  to  solve;  a with*  class  of  difficult 
and  important  problems,  although  much  more;  correlation  between  calculations  and 
experiments  is  needed. 


V.  THE  130-kbar  PHASE  CHANGE  IN  MARTENSITIC  STEEL 


In  recent  years,  a number  of  measurements  have  been  made;  of  the  high -pressure 
properties  of  iron  in  the  regions  of  the  a -►  e phase  transition.  These  indicate  that  at 
room  temperature  the  transition  occurs  in  the  region  of  130  khar  and  that  it  is  abarie, 
initiating  at  about  130  kbar  and  going  to  completion  at  pressures  > 170  khar.  The 
observed  specific  volume  change  associated  with  the  transition  is  A V = 0.0066  cm3/g 


^ D.  A.  Shockev,  “Ilyprrvdocity  Impact  and  Associated  Phenomena,"  SKI  Projccl  PYl  2151,  18  December  1973. 

^ “ A.  Slovens,  Personal  Communication,  SLA,  1973. 

1 3 

D.  A.  Shockey, "Hypervelocity  Impact  and  Associated  Phenomena,’'  SRI  ProjcclPYll  2151,  II!  December  197.’!. 


10 


and  the  compression  ratio  at  the  onset  of  the  transition  is  x = 0.943.  Here,  x is  the 
specific  volume,  v,  divided  by  the  specific  volume,  v , of  a - Fe  at  0°  K and  zero 
pressure.14 

The  phase  diagram  of  Fe  is  shown  in  Figure  5.  According  to  the  SLA  calcula- 
tions, the  shock  temperature  in  the  Fe  does  not  go  high  enough  to  drive  the  steel  into 
the  a phase.  The  shock  pressure  does  get  high  enough,  however,  to  drive  the  steel  into 
the  e phase  from  which  it  reverts,  on  cooling,  to  the  a phase,  with  a consequent  change 
in  grain  structure  and  mechanical  properties.  Direct  active  measurements  to  show  this 
change  have  not  been  made,  but  a combination  of  theory,  computations,  and  metallur- 
gical analysis  indicates  that  this  history  is  correct. 

Mosshauer  measurements15  show  that  the  magnetic  momentum  per  atom  fi  is 
practically  constant  for  the  a phase  up  to  the  pressure  of  its  disappearance  and  that 
the  e phase  is  non-ferromagnetic.  These  measurements  also  imply  that  the'  number 
N of  s-elcctrons  per  atom  is  invariant  with  pressure  for  the  a phase;  whereas,  for  the 
e phase,  Ns  decreases  with  increasing  pressure. 

Upon  the  reduction  of  pressure,  the  e phase  persists  down  to  a pressure  of  about 
45  kbar,  causing  hysteresis  in  the  compressional  behavior  of  iron  in  the  vicinity  of  the 
transition.  At  pressures  below  45  kbar,  the  e phase  disappears  completely.  Recent 
dynamic  measurements,  using  shock  waves,  have  been  made.16  These  confirm  the  non- 
ferromagnetic  nature  of  the  e phase  and,  hi  addition,  show  some  demagnetization  (an 
apparent  reduction  of  fi)  at  pressures  as  low  as  50  kbar. 

The  response  of  a metal  to  dynamic  loading  depends  upon  three  elements  which 
are  characteristic  of  the  material:  the  dynamic  constitutive  relation,  the  pressure- 

temperature  phase  relation,  and  the  dynamic  fracture  criteria.  In  the  absence  of  phase 
changes  or  fracture,  the  material  response  is  governed  completely  by  the  constitutive 
relation,  or  dynamic  equation  of  state.  The  constitutive  relation  is  characterized  as 
“elastic,”  “viscoelastic,”  “viscoelastic-plastic,”  etc.  It  depends  on  macroscopic  param- 
eters such  as  Youngs  modulus  and  the  viscosity  coefficient  which,  on  a microscopic 
scale,  arise  from  atomic  interactions  and  dislocation  processes. 

The  phase  diagram  determines  whether  a given  loading  history  will  bring  tin*  metal 
into  a pressure-temperature  regime  where  phase  transitions  are  possible.  If  this  occurs 
during  shock  loading,  multiple  shock  waves  may  result. 

t).  J.  Andrews,  “Equation  of  Slate  of  the  Alpha  and  Epsilon  Phase’s  of  Iron,”  Washington  State  University, 
WSU  SDL  70-05,  November  1970. 
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1).  J.  Pastine,  “A  Theory  of  the  a -*■  e transition  in  Fe  and  of  Possible  Higher  Pressure  Transitions  in  Fe  and  in 
tlie  Lighter  Elements  of  the  First  Transition  Series,”  NOL,  in  Metallurgical  Effects  at  High  Strain  Rates , Plenum 
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Figure  5.  Phase  diagram  of  iron. 


Figure  6.  Hugoniot  states  of  iron. 
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In  macroscopic  or  rheological  language*,  iron  can  he  described  as  a viscoelastic- 
plastic  material;  i.e.,  it  responds  elastically  until  yielding  (with  upper  and  lower  yield 
points)  and  is  strongly  rate -dependent.  The  constitutive  relation  thus  give's  the*  stress 
as  a function  of  strain  and  time  or,  alternatively,  can  he  written  as  a differential  equa- 
tion in  a,  e,  d,  and  e.  An  important  goal  of  present-day  research  is  to  develop  the* 
capability  to  derive  this  macroscopic  constitutive  relation  by  averaging  the*  effects  of 
microscopic  processes.  In  iron,  the*  density  of  mobile  el islocations  is  of  the*  order  of 
l()8/cin2,  which  is  large*  enough  to  make  a statistical  approach  valid. 

The*  elynainic  response!  of  iron  to  shock  loaeling  at  relatively  low  stre*ss  levels  is 
fairly  well  understood  on  the  basis  of  a dynamic  constitutive  re*l a tion  that  is  derived 
from  average  el isl e)cation  properties.  However,  it  is  not  yet  possible  uniquely  te)  derive* 
the  constitutive  relation  from  dislocation  properties  because  the*  ee>rrcct  dislocation 
model  based  on  statistical  averaging  of  dislocation  processes  has  not  been  completely 
determined.  Nevertheless,  it  is  to  be  expected  that  dynamic  experiments  selected  to 
distinguish  between  different  dislocation  models,  and  growing  understanding  of  the 
dislocation  mechanisms  themselves,  will  eventually  fill  in  most  of  the  remaining  gaps. 

As  noted  above,  a complicating  factor  in  the  dy  namic  response  of  iron  to  shock 
loading  is  the  presence  of  pressure-induced  phase  changes.  This  has  been  discussed 
briefly  before.  In  order  to  determine  whether  a given  loading  history  will  cause  a phase* 
change,  the  high-pressure  phase  diagram  of  iron  must  be  known.  The  pressure-induced 
phase  changes  in  iron  have*  many  formal  similarities  to  plastic  yielding.  Put  another 
way,  the  elastic-plastic  transition  with  strain  hardening  can  be*  regarded  as  a second- 
order  phase!  change,  since  it  involves  no  discontinuity  in  the  stress-strain  curve*  but  does 
involve  discontinuities  in  the  gradients.  The  kinetics  of  pi  as  tic  i ty  re!st  on  activation  and 
multiplication  of  dislocations,  and  the  kinetics  of  a crystallographic  phase!  transition 
are*  based  on  the  similar  hut  more  extensive  atennic  in e> tie>ns  which  change  the*  crystal 
structure  locally.  Thus,  the  constitutive  relations  can,  in  principle,  account  for  phase 
changes  in  a manner  analogous  to  the  way  they  handle*  plasticity. 

At  about  130  kbar,  a iron  begins  a transformation  to  a form  which  is  about  3 
percent  denser.  As  shown  in  Figure  5,  three  crystallographic  phases  of  iron  are  found 
to  be*  stable  in  different  regions  of  the  pressure-temperature  plane.  The  triple*  point  is 
at  about  1 15  kbar  and  775  K. 

The!  P-T  path  corresponding  to  high -velocity  impact  loaeling  of  iron  initially  at 
room  temperature  will  intersect  the  a - e coexistence  curve  at  a point  corresponding  to 
about  130  kbar.  The  effect  of  the  volume  compression  during  the  phase!  change  will 
cause  a shock  wave  of  initial  strength  between  130  and  330  kbar  to  split  into  two 
shock  waves  traveling  at  different  velocities.  Upon  reflection  from  free  surfaces,  the 
compressive  shock  waves  are  transformed  to  unloading  waves  which  may  intersect  to 
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produce  a tensile  stress  history  that  cause's  damage  and  failure.17  The  presence  of  the 
a ^ e phase;  change  can  result  in  extremely  steep  unloading  waves  or  “rarefaction 
shocks,’'  which  art*  effective  in  causing  fracture.  The  regime  in  the  P-T  plant*  reached 
during  loading  is  dependent,  of  course,  on  the  temperature  of  the  iron  prior  to  unload- 
ing and  on  whether  tin*  iron  is  stressed  and  heated  simultaneously  during  loading. 

VI.  INCIPIENT  SPALL  OR  FRACTURE  CRITERIA 

Dynamic  fracture  criteria  for  incipient  spall  have  been  considered  in  great  detail 
by  various  authors.18  This  is  discussed  tune  briefly.  Complete;  spall  and  back  face 
fragmentation  arc  discussed  in  Section  VII. 

In  the  ease  of  quasi-static  loading  of  brittle  materials,  it  is  often  possible  to  apply 
a simple;  fracture  criterion;  namely,  that  a critical  tensile  stress  exists  beyond  which 
failure  occurs.  Sample  geometry  and  the  effects  of  elastic  material  properties  and 
initial  crack  size  can  lx*  included  by  determining  the  fracture  toughness,  which  is 
related  to  tin*  maximum  stress  concentration  at  a crack  tip  allowable  before;  an  initially 
present  crack  propagates. 

This  simple  criterion  is  not  compatible  with  the  results  from  dynamic  experi- 
ments, however,  since  it  is  found  that  then;  is  a time  dependence  for  dynamic  failure; 
i.e*.,  tlx*  fracture  criterion  depends  on  the  entire  stress  history.  For  square  pulses,  this 
time  dependence  of  dynamic  fracture  has  been  shown  to  be  consistent  with  a relation- 
ship between  tensile  stress  and  tlx*  time  duration  of  tensile  stress,  At.  This  relationship 
is 

(o  - ow)x  At  = K 

where  a , X,  and  k an*  constants  specific  to  a particular  material.  If  X 1,  this  is  a 
simple  impulse*  criterion.  With  X = 2,  it  is  equivalent  to  an  energy  criterion.  With 
X = I,  the  equation  can  also  be;  related  to  other  time-dependent  fracture  criteria,  such 
as  stress-gradient  and  stress-loading  rati*,  which  arc  consistent  with  data  resulting  from 
explosive  tests.  These;  criteria  are  based  primarily  on  a simple  classification  of  damage 
into  categories  of  no  damage,  damage,  and  failure.  What  is  needed  is  a model  which 
relates  fracture  behavior  to  detailed  observations  of  the  damage  resulting  from  a given 
stress  history.  Such  correlations  an*  necessary  for  interpretation  of  tlx*  effects  of 
material  properties  on  dynamic  failure*  and  to  ensure  that  realistic  models  for  dynamic 
fracture*  are*  investigated.19 
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Another  approach  is  based  on  treating  damage  as  a new  phase  of  the  material, 
where  the  nueleation  and  growth  of  the  new  phase,  damage,  can  be  correlated  with 
macroscopic  stress  history  and  microscopic  material  properties.  In  other  words,  just 
as  phase  changes  can  be  treated  in  a manner  analogous  to  the  growth  of  plasticity,  the 
appearance  of  damage  can  also  be  treated  in  this  way.  A necessary  condition  is  that 
the  number  of  microcracks  be  sufficiently  largo  to  make  a statistical  approach  valid. 
The  steps  of  the  method  are:  (1)  shock  and  recover  specimens;  (2)  describe  the  size 
and  spatial  distribution  of  the  cracks  quantitatively;  (3)  achieve  experi mental  control 
so  that  the  cracks  can  be  stopped  in  different  stages  of  growth;  and  (4)  specify  the 
microscopic  stress  and  times-at-s tress  experienced  at  any  location  in  the  specimen.  If 
the  stress  history  is  sufficiently  uncomplicated,  then  it  is  possible  to  correlate  the 
observed  damage  to  the  stress  history  of  the  specimen  in  such  a way  that  nueleation 
and  growth  laws  of  general  validity  may  be  derived.  Once  these  factors  arc  known, 
quantitative  predictions  can  be  made  of  damage  caused  by  an  imposed  stress  history. 

This  approach  has  the  significant  advantage  of  bypassing  the  difficult  problem  of 
calculating  stress  distributions  around  microcracks.  Instead,  the  large  number  of 
mierocracks  present  allows  a statistical  approach,  and  the  crack  number  distribution 
functions  can  be  correlated  with  microscopic  stress  histories.  In  many  ways,  this  is 
analogous  to  the  derivation  of  constitutive  relations  from  average  dislocation  processes. 
The  large  number  of  dislocations  present  made  it  possible  to  use  a statistical  approach 
and  thus  avoid  the  problem  of  calculating  microscopic  stress  fields  around  individual 
dislocations. 


VII.  COMPLETE  ANI)  BACKPACK  SPALL 

As  discussed  in  Section  VI,  microcracks  and  voids  begin  to  form  when  the  tensile 
stress  reaches  the  “spall  threshold”  or  the  incipient  spall  strength.  For  practical  pur- 
poses, a certain  microcrack  size  or  density  is  assigned  to  this  value. 

As  the  tensile  stress  exceeds  the  spall  threshold,  the  microcracks  begin  to  coalesce 
until  a complete  spall  layer  is  formed.  For  the  impact  experiments  with  spherical 
projectiles  described  in  this  report,  the  surface  of  the  spall  layer  nearest  the  backface  is 
generally  parallel  to  the  backface  of  the  target  at  a distance  d from  the  backfaec.  For 
a givcm  type  of  projectile  d is  constant  over  the  velocity  range  from  3 to  7 km/s.  The 
spall  layer  has  a thickness  5 (perpendicular  to  the  backface),  which  increases  rapidly 
with  velocity  above  the  incipient  threshold.  In  addition  to  the  spall  layer  thickness, 
there  is  a backface  bulge  with  height  h,  which  is  equal  to  5.  As  long  as  the  spall  layer  is 
formed  near  the  backface,  h is  a useful  experimental  (or  empirical)  parameter  since  it 
can  be  measured  with  good  accuracy,  whereas  measurements  of  5 arc  sometimes 
ambiguous.  The  spall  layer  has  a diameter  D (parallel  to  the  backface)  that  increases 
slightly  with  impact  velocity.  It  grows  to  be  somewhat  greater  than  the  crater 
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diameter  (parallel  to  and  at  the  target  surface).  At  some  high  impact  velocity  the  baek- 
faee  “plugs,”  or  fragments;  i.e.,  material  of  thickness  d and  (approximate)  diameter  D 
is  ejected  from  the  baekfaee  either  as  a single  plug  or  as  a number  of  smaller  fragments. 
The  impaet  velocity  to  produee  incipient  spall  is  VQ  and  for  baekfaee  spall,  it  is  V^. 
As  will  be  seen,  appears  to  be  about  twice  VQ  for  the  experiments  described. 

Baekfaee  bulge  height  h is  plotted  on  linear  paper  as  a function  of  impaet  velocity 
for  three  different  eases  in  Figure  7.  It  is  seen  that  the  experimental  points  lie  rea- 
sonably well  on  smooth  eurves.  The  points  have  been  plotted  for  the  0.52-g  nylon 
sphere  on  the  1.25-cm  steel  target  on  semilog  paper  in  Figure  8.  The  points  fall  well  on 
a straight  line,  which  suggests  that  the  fractional  increase  in  h (or  5)  increases  with 
incremental  velocity ; i.e., 


A h = k 6 A V 

where  k is  a constant.  Upon  integration  this  becomes 

h = h ek<v-v°> 

O 

For  the  purposes  here,  ho  ean  be  assumed  to  correspond  to  the  bulge  height  at  inci- 
pient spall,  with  VQ  the  corresponding  impaet  velocity.  Empirical  values  for  k and  Vo 
are  given  in  the  following  table  for  the  three  eases  shown  in  Figure  7. 


Projeetile  Material 

Projeetile  Mass 
(g) 

Target  Thiekness 
(ern) 

h 

O 

(ern) 

k 

(s/km) 

vo 

(krn/s) 

Nylon 

0.52 

1.25 

0.01 

1.55 

2.77 

Nylon 

1.01 

2.5 

0.01 

1.40 

3.05 

112  O/Lexan 

0.35 

1.25 

0.01 

1.43 

3.23 

Admittedly,  while  much  more  experimentation,  analysis,  and  theory  are  needed, 
it  appears  that  the  above  results  may  provide  a good  working  criterion  for  complete 
spall  and,  possibly,  for  baekfaee  spall.  This  is  important  for  the  systems  engineer  who 
is  eoneerned  about  evaluating  certain  types  of  military  systems  at  minimum  eost.  If 
the  hypotheses  reported  here  are  true,  it  means  that  one  impaet  experiment  eliosen  for 
a given  target,  a given  projectile,  and  over  the  hyperveloeity  range  may  be  sufficient  to 
obtain  values  for  both  ineipient  spallation  and  baekfaee  spall. 
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PROJECTILE  TARGET 

C.  .52  GM  NYLON  1.25  CM  STEEL 


Figure  7.  Spall  damage  in  homogeneous  wrought  steel  plate. 
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Figure  8.  Backfaee  bulge  (or  spall  layer  thickness)  vs.  impact  velocity. 
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At  present,  2-d  computations  are  not  sufficiently  reliable  to  provide  the  neeessary 
predictions;  in  any  event,  they  are  prohibitively  expensive.  Furthermore,  a series  of 
impact  experiments,  such  as  those  shown  in  Figure  7,  is  also  expensive.  On  the  other 
hand,  one  shot  of  the  sort  described  is  relatively  inexpensive. 

In  order  to  justify  the  reported  hypotheses,  correlation  between  the  assumption 
that  h = 0.01  em  and  the  value  of  h corresponding  to  the  definitions  of  incipient  spall 
given  in  Section  VI  is  needed.  Based  on  experimental  results  and  analysis,  it  is  felt  that 
such  a practical  correlation  can  be  obtained. 

The  value  for  is  a different  matter.  Even  heuristic  theory  capable  of  pre- 
dicting is  not  available,  although  relevant  experiments  are  planned.  It  was  indi- 
cated previously  that  H20/Lcxan  and  nylon  spheres  are  somewhat  comparable  for  pro- 
ducing impaet-indueed  spallation.  It  is  seen  from  Figure  7 that,  for  a given  target  and 
for  plastie  and  liquid  spherical  projectiles,  V varies  inversely  as  the  projectile  mass. 
The  relation  between  target  thickness  and  projectile  eharaeteristies  is  not  elear, 
although  some  crude  empirical  formulas  may  be  possible. 

VIII.  SUMMARY  ANI)  CONCLUSIONS 

The  results  of  extensive  hypervclocity  impact  experimentation  on  martensitie 
steel  armor  have  been  summarized.  It  has  been  suggested  that  certain  types  of  austeni- 
tic steel  may  be  less  subject  to  spallation  damage  than  martensitic  steel.  Damage  cri- 
teria from  complete  and  baekface  spall  have  been  discussed,  and  heuristic  or  empirical 
theory  for  predicting  baekfaee  spall  in  martensitie  steel  armor  has  been  suggested.  The 
main  conclusion  of  this  work  is  that  the  a - e phase  change  for  martensitie  steel  is 
important  in  predicting  or  interpreting  projectile  impact  damage.  Other  conclusions 
are  more  general  and  possibly  not  applicable  for  present  purposes.  For  example,  a 
preliminary  conclusion  is  that  hypervclocity -impact-induced  baekfaee  spall  is  a viable 
damage  mechanism.  Another  preliminary  conclusion  is  that  efficient  spallators  capable 
of  optimizing  baekfaee  spall  arc  feasible.  The  main  conclusion  is  that  more  experi- 
ment, theory,  analysis,  and  computations  are  needed. 
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